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ABSTRACT. Published infrageneric classifications of Coccothrinax are based only on Cuban species. To 
provide an inclusive categorization of provisional infrageneric groupings for easy recognition in future field 
and herbarium studies, we compared all species in the genus with statistical, phenetic, and cladistic tech- 
niques. Statistical analyses indicated that 22 of 160 characters extracted from the published literature would 
be useful for phenetic and cladistic analyses. Phenetic analyses produced three distinct clusters of taxa. 
Cladistic analyses suggested that two of the phenetic clusters are paraphyletic and the third is monophyletic. 
The two paraphyletic clusters were also more plesiomorphic than the monophyletic cluster. Cladistic 
analyses further indicated that most of the 22 characters were highly homoplasious, and only the following 
were reliable for classificatory purposes at the infrageneric level: loss of transverse veins in the leaf, elongation 
of the leaf sheath strands, fusion of sheath strands into spine-like processes, and tan, corky fruits. The 
recognizable groups are as follows: Argentea Group (1 complex, 6 spp., paraphyletic); Argentata Group (2 
complexes, 10 spp., probably paraphyletic) — C. argentata complex (8 spp., paraphyletic), C. crinita complex 
(2 spp., monophyletic); and Pauciramosa Group (4 complexes, 31 spp., probably monophyletic)—C. mira- . 
guama complex (11 spp., probably monophyletic), C. pauciramosa complex (11 spp., paraphyletic), Haitiella 


complex (2 spp., monophyletic), species incertae sedis (7 spp.). 


As a result of recent name changes and an 
absence of detailed comments on species rela- 
tionships, the taxonomy of Coccothrinax Sargent 
has become complex and poorly resolved. The 
number of recognized taxa in the genus is unsta- 
ble and relatively large. The total number of spe- 
cies varies from 20 cited by Moore (1973) to 37 
listed by Glassman (1972). Since 1966, numer- 
ous name changes have been proposed and 14 
new taxa have been described (Borhidi & Muniz, 
1971, 1972, 1985; Borhidi et al., 1978; Muniz 
& Borhidi, 1982a, 1982b; Quero, 1980; Read, 
1966a, 1966b, 1980). Additionally, the only in- 
frageneric classification (Leon, 1939) was based 
exclusively on Cuban taxa. Leon’s (1939) scheme 
has been variously modified by Muniz and Bor- 
hidi (1982b). These authors recognize two sec- 
tions and four subsections, but do not elaborate 
on relationships among the groups, except as im- 
plied by the taxonomic rank of the groups. Spe- 
cies relationships within the groups are also not 
given, except as isolated comments noting affin- 
ities among two or three species. Muniz and Bor- 
hidi’s modifications of León's scheme are also 
based on an analysis of only the Cuban taxa. 


* Present address: General Biology Program, 301 
Hesler Biology Bldg., University of Tennessee, Knox- 
ville, Tennessee 37996-1110. 


The purpose of this study is to provide a hier- 
archial context for future detailed populational 
and species-level systematic studies, and to pro- 
vide insight into the relationships within the ge- 
nus. We examined the available literature data, 
augmented with field, herbarium, and garden 
studies, and developed an inclusive preliminary 
scheme that is hierarchially more complex than 
that of León ex Muniz and Borhidi (1982b). 


MATERIALS AND METHODS 


Data set. Literature data, including previously 
published keys, descriptions, photographs, and 
drawings, were used to characterize the taxa for 
comparison. A list of recognized taxa is given in 
TABLE 1 and a list of literature consulted is cited 
in Nauman and Sanders (1991). The literature 
data were augmented with single or multiple 
specimen samples from 44 authenticated living 
accessions at Fairchild Tropical Garden (FTG), 
the herbarium, and field studies by one or both 
authors in Florida, Hispaniola, Mexico, and 
Trinidad (specimens cited in Nauman & Sanders 
(1991), APPENDICES I and II). 

Character analysis. More than 160 characters 
used by previous authors to distinguish taxa and 
additional characters considered to be of poten- 
tial value for indicating systematic relationships 
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TABLE 1. 


Cn RB WN = 


SELBYANA 


List of recognized taxa. Numbers refer to 
taxa in FIGURES 2 and 3, and APPENDIX II. Lit- 
erature sources and synonymy are given in Nau- 
man and Sanders (1990). Infraspecific taxa have 
been omitted at this time. 


. C. acunana León 

. C. alexandri León 

. C. alta (Cook) Beccari 

. C. argentata (Jacquin) Bailey 

. C. argentea (Lodd. ex J. A. & J. H. Schultes) Sar- 


gent ex K. Schum. 


. C. baracoensis Borhidi & Mufiiz in Mufiiz & Bor- 


hidi 


. C. barbadensis (Lodd. ex Martius) Beccari 

. C. bermudezii León 

. C. borhidiana Muñiz in Borhidi et al. 

. C. camagueyana Borhidi & Mufiiz in Mufiiz & 


Borhidi 


. C. clarensis León 
. C. concolor Burret 
. C. crinita (Griseb. & Wendl. ex Wright in Sauvalle) 


Beccari 


. C. cupularis (León) Borhidi & Muñiz 

. C. ekmanii Burret 

. C. elegans Muñiz & Borhidi 

. C. fragrans Burret 

. C. garciana León 

. C. gracilis Burret 

. C. guantanamensis (León) Mufiiz & Borhidi 
. C. gundlachii León 

. C. hioramii León 

. C. inaguensis R. W. Read 

. C. jamaicensis R. W. Read 

. C. leonis Muñiz & Borhidi 

. C. litoralis León 

. C. macroglossa (León) Mufiiz & Borhidi 

. C. microphylla Borhidi & Muñiz in Muñiz & Bor- 


hidi 


. C. miraguama (H.B.K.) León 
. C. moaensis (Borhidi & Muñiz) Muñiz in Muñiz 


& Borhidi 


. C. montana Burret 

. C. munizii Borhidi in Borhidi & Muñiz 

. C. muricata León 

. C. nipensis Borhidi & Muñiz in Muñiz & Borhidi 
. C. orientalis (León) Mufiiz & Borhidi 

. C. pauciramosa Burret 

. C. proctorii R. W. Read 

. C. pseudorigida León 

. C. readii Quero 

. C. rigida (Griseb. & Wendl. ex Griseb.) Beccari 

. C. salvatoris León 

. C. savannarum (León) Borhidi & Muñiz in Muñiz 


& Borhidi 


. C. saxicola León 

. C. spissa Bailey 

. C. victorinii León 

. C. yunquensis Borhidi & Muñiz in Muñiz & Bor- 


hidi 


. C. yuraguana (A. Rich.) León 
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were surveyed for taxonomic utility. Patterns of 


character variation were analyzed, and coding 


procedures were determined from studies of four 
disjunct Floridian populations of C. argentata 
and living accessions of C. argentea, C. argen- 
tata, C. miraguama, C. readii, C. yuraguana, and 
C. proctorii at FTG (unpubl., data and results 
available upon request from C.E.N.). Herbarium 
material at FTG was used to supplement the stud- 
ies of live material. 

R-mode analyses, including multiple discrim- 
inant analysis, principal component analysis, and 
other statistical procedures, were applied to par- 
ticular characters and garden populations to as- 
sess their usefulness for intertaxon discrimina- 
tion. Several rejection criteria were used to 
evaluate each character. 1) The character had to 
exhibit some variation within or among the OTU's 
(operational taxonomic units, here considered 
equivalent to a named taxon, TABLE 1). This pro- 
cedure emphasizes character differences among 
the oTu’s. Additionally, where intra-OTU varia- 
tion was greater than inter-OTU variation, i.e., 
where F-1 weighting (sensu Adams, 1975) would 
give zero weight to the character, differences 
among OTU’s would not be detectable; these char- 
acters were eliminated. 2) The character had to 
be measurable on a large proportion of the avail- 
able specimens. 3) The character had to exhibit 
measurement error smaller than the differences 
in that character among any two or more OTU's. 
Using these criteria and combining logically or - 
empirically correlated characters into a character 
complex (e.g., tan corky fruits), the initial 160 
characters were narrowed to 22 (APPENDIX I). 

Character state coding for phenetic analyses 
was determined with statistical summaries, in- 
cluding a series of basic descriptive statistics. For 
most characters, modal states were considered 
representative due to their higher predictive val- 
ue (cf. Nauman, 1982). This was especially im- 
portant for characters exhibiting non-normal 
distributions. A character state distribution was 
considered non-normal if the third or fourth cen- 
tral moment statistics indicated significant skew- 
ness or kurtosis, P « 0.05. Modal states for con- 
tinuous characters were obtained by partitioning 
the range into ten class intervals. In some in- 
stances, continuous characters were broken into 
classes delimited by natural breaks in the char- 
acter state frequency distributions using a meth- 
od similar to that of Almeida and Bisby (1984). 
Some characters were represented by maximal 
values for each taxon since these values were 
shown from analysis of field or garden popula- 
tions to be discriminating only at the upper end 
of their ranges. Additionally, multivariate anal- 
ysis indicated that the upper range limits for cer- 
tain characters represent the only comparable 
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developmental states. Basic statistical summa- 
ries were carried out with programs written by 
one of us (C.E.N.) on an IBM PC. The raw data 
matrix is given in APPENDIX II. 

Phenetic analyses. The phenetic analyses con- 
sisted of: 1) generating the general similarity co- 
efficient of Gower (1971) for all pairwise com- 
binations of OTU's; 2) cluster analysis (CA) of the 
similarity matrix with four pair-group clustering 
methods, UPGMA—an unweighted method with 
arithmetic averages, WPGMA—a weighted meth- 
od with arithmetic averages, weighted centroid, 
and unweighted centroid methods; and 3) Prin- 
cipal Coordinate Analysis (PCORD) on the simi- 
larity matrix with the method of Gower (1966). 
The WPGMA and centroid clustering were carried 
out with the Multivariate Statistical Package (W. 
L. Kovach, Indiana Univ.). Computation of the 
similarity matrices, UPGMA, and PCORD were car- 
ried out with a phenetic program package written 
by one of us (C.E.N.) on an IBM PC. 

Cladistic analyses. For cladistic analyses, poten- 
tial homology was estimated by similarities in 
positional and developmental relationships. Since 
such similarities can arise independently, con- 
gruence with the entire data set was considered 
the stronger test of homology. Polarities were 
assessed largely with the outgroup substitution 
approach of Donoghue and Cantino (1984). In 
most cases, Thrinax and Zombia were consid- 
ered the immediate sister groups of Coccothrinax 
(cf. Uhl et al., 1990). When polarity assessment 
was conflicting, other more distantly related cor- 
yphoid genera were utilized. Where outgroup 
substitution was inappropriate for resolving po- 
larities (e.g., when plesiomorphic states were am- 
biguous or non-existent in the outgroup) we used 
the “ontogenetic” criterion and compared seed- 
lings of a number of taxa grown at FTG. Addi- 
tional details of the polarity assessments are giv- 
en in APPENDIX I. Computations were carried out 
on an IBM PC with PAUP (version 2.4.1) obtained 
from D. L. Swofford (Illinois Natural History 
Survey). 


RESULTS 


Phenetic analyses. The results of both CA (FIGURE 
1) and PCORD (FIGURE 2) suggest three main phe- 
netic clusters, here designated as the Argentata, 
Argentea, and Pauciramosa groups. Only the 
phenogram resulting from UPGMA clustering of 
47 OTU's is presented (FIGURE 1) because the other 
three CA procedures were no more informative. 
In PCORD (FIGURE 2), positive scores on Axis 1 
are associated with smaller leaves (segment length 
and number, and palman length, 1.e., segment 
connation) and more open spacing of sheath 
strands. Positive scores on Axis 2 are associated 
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with a lack of transverse veins, thicker sheath 
strands, and free strand tips formed by strand 
fusion (spine-like sheath tips). Positive scores on 
Axis 3 are associated with slender caudices. 

The major discrepancy between the results of 

CA and PCORD is the placement of C. acunana, 
C. gundlachii, and C. montana. In CA, they arise 
at low levels of similarity, being clearly associ- 
ated with none of the more distant clusters. In 
PCORD, they appear to be intermediate between 
the Argentata and Pauciramosa groups. These 
species recombine character states from these two 
groupings and the conflicting placement reflects 
how these two algorithms portray inter-OTU and 
inter-cluster relationships; CA tends to portray 
inter-OTU relationships more accurately, and 
PCORD tends to portray inter-cluster relationships 
more accurately. 
Cladistic analyses. Approximately 100 minimal 
length trees were obtained with PAUP using the 
MULPARS option with global branch swapping and 
rooting by the Lundberg method. Since polarities 
were previously assigned by outgroup compari- 
sons, we scored plesiomorphic states as 0 and 
apomorphic states as 1 or 2. Thus, the Lundberg 
rooting simulated the presence of a hypothetical 
ancestor possessing the plesiomorphic condition 
for all characters. The trees have a branch length 
of 69 and a Consistency Index of 0.41. A con- 
sensus tree of Adams for those 100 trees is given 
in FIGURE 3; it has a branch length of 70 and a 
Consistency Index of 0.40. 

The character conflict in C. acunana, C. gund- 
lachii, and C. montana is handled differently by 
cladistic procedures. Most of the conflicts in C. 
acunana and C. montana are in plesiomorphic 
character states; hence, the taxa are placed at the 
base of the tree. Those of C. gundlachii are pri- 
marily apomorphic in the context of its parsi- 
moniously placed neighbors; hence, it is placed 
at the top of the tree. 

The Argentea phenetic group remains coher- 
ent but is a paraphyletic assemblage whose com- 
mon characters are plesiomorphic. The Argen- 
tata phenetic group similarly remains coherent 
and is also paraphyletic but intermediate in po- 
sition, and no apomorphic character states unite 
the orv's. The Pauciramosa phenetic group is an 
advanced monophyletic clade defined by char- 
acters of few segments and short palmans, both 
of which are homoplasious. 

Only ten apomorphous character changes (2- 
1, 4-2, 6-2, 7-1, 8-1, 15-1, 17-1, 18-1 [-2], 20-1, 
and 22-1) of the possible 28 do not exhibit ho- 
moplasy. All of these except 6-2 and 17-1 are 
autapomorphies for species or closely related 
species pairs. Character state 6-2 (long sheath 
spines) unites nine taxa; 17-1 (no transverse veins) 
unites 39 taxa. Therefore, almost all of the hier- 
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FIGURE 1. 


Phenogram resulting from UPGMA cluster analysis of 47 oTU’s. 
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FIGURE 2. Three dimensional diagram of the results of Principal Coordinate Analysis of 47 oru's (OTU 


numbers are given in TABLE 1). 


archial structure of the tree is based on homopla- 
sious characters with individual consistency val- 
ues of 0.5 or less. If the tree structure is reduced 
to the two consistent apomorphies, we find 1) a 
clade of species lacking transverse veins arising 
from a primitive complex, and 2) within that 
clade, a secondary clade with long sheath spines. 


DISCUSSION AND CONCLUSIONS 


Limitations of the data. As one might expect, the 
data that can be obtained from species descrip- 
tions have limitations for broad comparisons and 
for numerical analyses. If the species are known 
only or primarily from type material, descrip- 
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FIGURE 3. Consensus tree of Adams for 100 trees resulting from parsimony analysis of 47 oTU’s. Numbers 
represent character state changes at each position, where the first number refers to the character and the second 
number refers to the resulting character state (see APPENDIX I). 


tions tend to be incomplete. Thus, many of the 
160 characters initially examined were not com- 
parable across the genus. This is especially true 
for reproductive characters. (Note that only four 
reproductive characters are included in APPENDIX 
I and only one of these concerns the flower.) 
Vegetative characters are prone to parallelisms. 

Most of the comparable vegetative characters 
from the literature are described quantitatively. 
When an attempt is made to convert the reported 
measures into discrete states, however, often only 
the extremes are distinguishable as morpholog- 
ical gaps. Hence, in our study, eight of the ten 
non-homoplasious characters were unique to a 
single species or species pair. As a result, the 
overall phenetic structure was derived from the 
quantitative vegetative characters. This is not 
problematic for most phenetic algorithms, but it 
does mean that the phenetic groups derived from 
such characters may not accurately reflect evo- 
lutionary relationships. In a cladistic context, 
these same characters exhibit relatively high de- 
grees of homoplasy. P. F. Stevens (unpubl.) ar- 


gues that the binary or ordinal coding of quan- 
titative characters may not be valid unless the 
intraspecific variation is well documented and 
there are absolute gaps between the ranges of 
species scored for different character states. Be- 
cause descriptions of new species often do not 
fulfill the first criterion, these types of data may 
be further limited for cladistic analysis. 

We conclude that the clusters and hierarchy 
obtained from these data should not be directly 
taken at face value. Rather, they should be used 
as guides for direct observations of known as well 
as previously under-utilized characters. Further- 
more, characters that were eliminated from con- 
sideration by one of the primary rejection criteria 
in the context of the entire genus take on a new 
significance in the context of individual groups 
resulting from the phenetic and cladistic analy- 
ses. Thus, we were led to re-evaluate the char- 
acters of a wider sample of specimens for ab- 
solute gaps rather than statistical differences. For 
example, the presence of transverse veins (char- 
acter 17) in the Argentea phenetic grouping dif- 
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TABLE 2. Classification of Mufiiz and Borhidi (1982b). Taxa indicated by numbers in TABLE 1. Asterisks 
indicate taxa outside Cuba and not treated by them. 








Coccothrinax Sargent 


Section Coccothrinax 
Subsection Coccothrinax: 


1 2,3" 4*..5*. 7*5; 8, 14, 16, 17; 19*. 20. 22, 29*. 24*. 


25-30, 31*, 33, 35, 37, 39, 41, 43, 44*, 45-47. 


Subsection Haitiella (L. H. Bailey) Muñiz: I5*. 32. 
Section Longispadiceae León ex Muñiz 
Subsection Multiramosae León ex Muñiz: 10, 13, 21. 


Subsection Pauciramosae León ex Muñiz: 


6, 9, 11, 12*, 18, 34, 36, 38, 40, 42. 





fers dramatically from the absence of them in 
the other two phenetic groupings. (A minor ex- 
ception is noted in Nauman and Sanders (1991) 
which may be due in part to garden hybridiza- 
tion.) The use of this character as a cladistic char- 
acter is substantiated. On the other hand, direct 
observations of homoplastic or quantitatively 
varying characters show that personal interpre- 
tation more strongly influences the way previous 
authors described them or different workers score 
them. In this context, a comparison of our results 
with the published classifications provides ad- 
ditional insights. 

Published classifications. The scheme of classi- 
fication for the genus Coccothrinax published by 
Muniz and Borhidi (1982b) (TABLE 2) recognizes 
two sections, section Coccothrinax and section 
Longispadiceae. Section Coccothrinax is char- 
acterized by short, cernuous to pendent inflores- 
cences before fruit maturity, 1 to 6 primary in- 
florescence branches, and partly membranous 
external bracts. As circumscribed by Muniz and 
Borhidi (1982b), the section contains two sub- 
sections, subsection Coccothrinax and subsec- 
tion Haitiella. Subsection Haitiella (previously 
Haitiella Bailey) is distinguished by its densely 
muricate ovaries and fruits, by leaves that are 
less than !^ orbicular, by short (less than 4 dm 
long) segments, and by short (less than 2—2.5 dm 
long) petioles. These character states occur scat- 
tered in other subsections but are typical of mem- 
bers of subsection Haitiella. Subsection Coc- 
cothrinax is characterized by smooth ovaries and 
fruits. 

Muniz and Borhidi (1982b) characterized sec- 
tion Longispadiceae by elongated and ascending 
inflorescences with rigid external bracts and di- 
vided the section into two subsections, subsec- 
tion Pauciramosae and subsection Multiramo- 
sae. Subsection Pauciramosae is characterized 
by small- to medium-sized plants with 1 to 7 
primary inflorescence branches, and subsection 
Multiramosae is characterized by robust plants 
with 8 to 10 primary inflorescence branches. 
Comparison of schemes. Section Longispadiceae 


corresponds to the Pauciramosa phenetic group, 
and section Coccothrinax corresponds to the Ar- 
gentea and Argentata groups combined. The ma- 
jor exceptions are: 1) Muniz and Borhidi (1982b) 
placed C. borhidiana, C. crinita, and C. cama- 
gueyana in section Longispadiceae because they 
possess long, horizontal to ascending inflores- 
cences with rigid bracts; 2) because C. ekmanii 
and C. munizii have short, cernuous inflores- 
cences with partly membranous bracts, they are 
placed in section Coccothrinax as subsection 
Haitiella; and 3) several species they excluded 
from subsection Pauciramosae fell in the Pau- 
ciramosa phenetic group. The need to segregate 
species of the Argentea group from section Coc- 
cothrinax would not have been obvious to Cuban 
authors because only two of the taxa occur in 
Cuba (TABLE 2). 

Prior to our numerical analyses, we were un- 
able to confirm the distinctions between rigid and 
partly membranous bracts and between ascend- 
ing and cernuous inflorescences; these characters 
were therefore excluded from the numerical anal- 
yses. Although Muniz and Borhidi (1982b) relied 
heavily on inflorescence characters, their scheme 
is remarkably similar to our results. We believe 
that in addition to their use of the inflorescence, 
they were intuitively influenced by the same 
characters used in our analyses. Indeed, when 
subsection Haitiella is placed within subsection 
Pauciramosae, the latter can be defined by free 
sheath strand tips (sheath spines) over 1 cm long 
(character state 6-1) as in FIGURE 3. 

The reliance on inflorescence posture and 
length, however, appears to have created artifi- 
cial associations. For example, Mufiiz and Bor- 
hidi (1982b) recognized that C. borhidiana and 
C. crinita share the same type of unusual free 
strand tips. However, because the former has 
shorter inflorescences they were forced to place 
the two in different subsections, despite their ob- 
vious relatedness in other characters that are more 
conservative than inflorescence length. Further- 
more, based on observations of seedlings, we be- 
lieve that the free strand tips in C. crinita are 
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TABLE3. Informal classification of Coccothrinax, de- 
duced from the numerical and cladistic analyses 
(see TABLE 1). “+” indicates “more or less." 


1. ARGENTEA GROUP: Most primitive, probably 
paraphyletic. Cuba to Trinidad, only 1 complex. A. 
A. ARGENTEA COMPLEX: Stems ca. 10-30 cm 

or more in diameter; leaf sheath strands fine, 
hairlike, usually much less than 0.5 mm thick, 
tightly woven; leaf blade shallowly plicate in one 
plane, often with segments distally pendent; seg- 
ments usually tapering gradually to tip; longi- 
tudinal veins widely spaced, interconnected by 
distinct transverse veins; inflorescence bracts + 
tubular. Taxa: 3, 5, 7, 20, 22, 44. 


2. ARGENTATA GROUP: Intermediate advance- 
ment, probably paraphyletic. Stem mostly 10-20 
cm in diameter; leaf sheath strands fine, usually 
hairlike, usually less than 0.5 mm thick, usually 
tightly woven; leaf blade shallowly plicate, often 
broadly undulating or funnel shaped, often with seg- 
ments distally pendent; segments usually tapering 
gradually to the tip; longitudinal veins numerous, 
closely spaced, not interconnected by distinct trans- 
verse veins; inflorescence bracts + tubular. North- 
western Caribbean, 2 complexes. 

A. ARGENTATA COMPLEX: Probably paraphy- 
letic. No additional defining characters. Taxa: 
4, 17, 23, 24, 26, 37, 39, 45. 

B. CRINITA COMPLEX: Probably monophylet- 
ic. Free strand tips of leaf sheath due to elon- 
gation. Taxa: 9, 13. 


3. PAUCIRAMOSA GROUP: Advanced, probably 
monophyletic. Stem usually 3-10 cm in diameter; 
leaf sheath strands coarse, usually woody, 1-4 mm 
or more thick (wiry and only 0.5 mm thick in a few 
taxa of the Miraguama complex); leaf blade deeply 
plicate in one plane, usually with segments distally 
held stiffly straight (hence, noticeably pinwheel 
shaped); segments usually abruptly contracted and 
producing shoulder-shaped acuminations distally; 
longitudinal veins numerous, closely spaced, not 
interconnected by distinct transverse veins; inflo- 
rescence bracts tubular to dilated. Cuba, with some 
taxa in Hispaniola, 3 or more complexes. 

A. MIRAGUAMA COMPLEX: Monophyletic or 
paraphyletic. Leaf sheath strands in 3 layers (or 
occasional apparent reversal to 2 layers), + 
tightly woven, free spiniform strand tips poorly 
developed to ca. 1 cm long or lacking; leaf seg- 
ments variable in number, length, and conna- 
tion. Taxa: 2, 6, 8, 10(?), 14, 25, 27, 29, 30, 35, 
47. 

B. PAUCIRAMOSA COMPLEX: Probably para- 
phyletic. Leaf sheath strands in 2 layers, + 
loosely woven, free spiniform strand tips well 
developed (over 2 cm long); leaf segments usu- 
ally fewer than 35, short, and weakly connate. 
Taxa: 11, 18, 21, 28, 33, 34, 36, 38, 41, 42, 46. 

C. HAITIELLA COMPLEX: Monophyletic. Leaf 
sheath strands in 2 layers, + tightly woven, free 
spiniform strand tips well-developed; leaves cu- 
neate with segments fewer than 30, short, weak- 
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TABLE 3. Continued. 


ly connate; ovaries muricate, fruit tan, corky 
(autapomorphies in Coccothrinax). Taxa: 15, 32. 

D. MISCELLANEOUS TAXA: Leaf sheath 
strands in 2 layers, + loosely woven, free spi- 
niform strand tips present or not; leaf segments 
variously but typically fewer than 35, short, and 
weakly connate. Taxa: 1, 12, 16, 31, 40, 43. 
(Included here is C. scoparia, once treated as a 
synonym of C. argentea by Bailey (1939). It was 
excluded when we conducted the numerical and 
cladistic analyses, but later study of authentic 
material showed it to be distinct and most re- 
lated to C. saxicola.) 


homologous with the fused strand tips in the 
Argentata group and not with the free spinelike 
strand tips in subsection Pauciramosae. 

Of special interest are the taxa positioned in 
subsection Coccothrinax by Muniz and Borhidi 
(1982b) and which clustered in the Pauciramosa 
group in our analyses. These taxa have sheath 
strands consistently more than 0.5 mm thick (wiry 
to woody in consistency), mostly loosely woven 
sheath strands, and leaf segments that are short, 
few per leaf, and only shortly connate. Also, the 
taxa centered around C. yuraguana share the 
apomorphy of three sheath layers. However, these 
same characters are variously recombined in the 
taxa centered around C. miraguama (C. alex- 
andri, C. cupularis, C. leonis, and C. macroglos- 
sa). In contrast, these taxa are aligned with the 
Argentata phenetic grouping by both Muniz and 
Borhidi (1982b) and our analyses. 

With respect to these same characters, C. mi- 
raguama is sufficiently variable to be partitioned 
into subspecies. Coccothrinax miraguama subsp. 
arenicola approaches C. yuraguana, a species with 
wiry sheath strands in three layers, but differs 
from the latter in more numerous leaf segments 
that are more strongly connate. Subspecies hava- 
nensis and roseocarpa have thick, woody sheath 
strands in three layers as in C. orientalis but also 
have larger, more numerous, and more connate 
segments; subsp. roseocarpa is the more transi- 
tional of the two subspecies. Subspecies mira- 
guama is nearly identical to subsp. havanensis, 
but has the sheath strands typically in two layers. 
However, in this subspecies the distinction be- 
tween two and three sheath strand layers breaks 
down; three layers are often present proximally 
and two layers distally on the sheaths. Short, free 
strand tips about 1 cm long are usually present 
in C. cupularis, C. macroglossa, and some plants 
of C. miraguama, which suggests genetic simi- 
larity to other taxa in the Pauciramosa phenetic 
group. Thus, it appears that C. miraguama and 
associated taxa would be better placed in the 
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Pauciramosa group near C. yuraguana. Indeed, 
R. W. Read (pers. comm.) has opined that the 
two taxa may be conspecific. Such a move, in 
our opinion, makes the Argentata phenetic group 
more homogeneous and places in the Paucira- 
mosa phenetic group all the taxa with similar 
character state combinations. This change would 
also reduce the homoplasy in certain characters, 
e.g., the number of sheath layers (character 3) 
and the presence of spinelike free sheath tips 
(character 6). The homoplasies that would re- 
main, or become more homoplasious as a result 
of this action, are quantitative vegetative char- 
acters (leaf segment number (character 11), length 
(character 12), and fusion (character 13)). 
Thus, it is clear that reliance on published data 
provides only a few significant classificatory 
characters, as well as numerous characters that 
are weak indicators of relationships. In our anal- 
yses, the latter tend to swamp out the former. 
However, by completing multiple analyses and 
by comparing the results with the published in- 
tuitive classifications and with living plants, we 
are able to draw two significant conclusions. First, 
the distinction between more reliable and less 
reliable characters is clarified. Second, additional 
morphological features that suggested patterns of 
affinities can be correlated with the reliable char- 
acters to produce eclectic, easily recognizable 
groups of intermediate size (TABLE 3). These 
eclectic groups are temporarily useful in assisting 
workers to focus efforts in field studies and on 
determining species limits. Future studies should 
concentrate on character data, both vegetative 
and reproductive, that can be compared among 
all species before meaningful hypothesis testing 
can begin. Our arrangement of taxa in TABLE 3 
provides a framework for such studies as well as 
preliminary but readily testable hypotheses. 
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SELBYANA 


APPENDIX I. CHARACTERS AND CHARACTER STATE 


Copinc. All character states coded as 0 are con- 
sidered plesiomorphic; apomorphies are coded 1 
or 2. The method of polarity assessment is given 
in parentheses following the character states. 


. Caudex maximum diameter (excluding persistent 
leaf sheaths and petiole bases). O — generally «2 
dm; 1 = generally 22 dm. (outgroup). 

. Caudex up to 3 dm and ventricose. 0 = no; 1 = 
yes. (outgroup). 

. Sheath layers. The leaf sheath is composed of a 
network of enlarged fibrovascular bundles. These 
are the “sheath fibers" of other students of Coc- 
cothrinax. We here follow Tomlinson (1964) in 
calling them strands. The strands are arranged in 
layers where those within a layer are parallel to 
one another and more or less perpendicular to 
those of adjacent layers. Strands within a layer are 
held together by a reticulum of decidedly smaller 
cross-strands that branch at right angles to, and 
interconnect with, adjacent main strands. The 
number of layers was determined from mature 
sheaths near or at the apex of the sheath body, or 
base of the free tip if present. 0 = 2 or 2 with a 
3rd indistinct layer; 1 — 3 layers. (outgroup). 

. Leaf sheath strand thickness, as the average width 
of the thicker strands in the body of the leaf sheath. 
0 = up to 2 mm wide; 1 = 2.1-4.9 mm wide; 2 = 
5 or more mm wide. (outgroup). 

. Leaf sheath strand spacing at the sheath body. O 
= tight (i.e, <1 x the strand thickness between 
strands); 1 = open. (outgroup). 

. Length of free leaf sheath strand tips (those due to 
fusion of bundles and lack of cross connections) 
at the sheath tip. 0 = «1 cm long; 1 = 1-2.9 cm 
long; 2 = =3 cm long. (ontogeny). 

. Length of free leaf sheath strands (those due to 
elongation of the sheath tip and fragmentation of 
the sheath, i.e., splitting of the strands from each 
other) greater than 10 cm. 0 = no; 1 = yes. (on- 
togeny). 

. Free leaf sheath strands with uncinate tips. 0 = 
no; l = yes. (ontogeny). 
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Continued. 


. Petiole length routinely less than 2 dm long. 0 — 


no; 1 = yes. (outgroup). 


. Lamina outline less than 2 orbicular. 0 = no; 1 = 


yes. (outgroup). 


. Segment number. 0 = 55 or more; 1 = 35 to 54; 


2 = 34 or less. (outgroup). 


. Segment length. 0 = 4.1-6.9 dm; 1 = =7 dm; 2 = 


<4 dm. (outgroup). 


. Greater palman length. The palman is that portion 


of the leaf consisting of fused segments and ex- 
cludes the free segment apices. Greater palman 
was measured as the length of the fused region 
between the central pair of segments. 0 = 21.3 
dm; 1 = <1.3 dm. (outgroup). 

Segments widest past the sinus. 0 = yes; 1 = in- 
termediate or about at the sinus; 2 = no. (out- 
group). 

Lorica remnants generally present. 0 = no; 1 = 
yes. (outgroup). 

Abaxial lamina indumentum color. 0 = concol- 
orous; 1 = grayish, white, or silvery. (ontogeny). 
Transverse veins obviously present. These are veins 
running perpendicular to the long axis of the seg- 
ment. Such veins were considered present only if 
a well-developed reticulum is exhibited, 1.e., one 
that can be easily seen with the unaided eye or a 
handlens. 0 = yes; 1 = no. (outgroup). 

Adaxial hastula (ligule) apex bifid. 0 = no; 1 = 
sometimes; 2 = characteristically. (outgroup and 
ontogeny). 

Rachilla length. 0 = «24 cm; 1 = =24 cm. (out- 
group). 

Anther length. 0 = <5 mm long; 1 = =5 mm long. 
(outgroup). 

Ovary surface muricate. 0 = no; 1 = yes. (ontog- 
eny). 

Exocarp tan at maturity. 0 = no; 1 = yes. (outgroup 
and ontogeny). 
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OTU numbers are those listed in 
12 13 14 15 16 17 18 


l 
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0 0 0 0 0 0 0 0 0 0 


l 


unknown or unavailable, the most parsimonious states were inserted by PAUP; these are indicated by 


TABLE 1, character numbers are those listed in APPENDIX I. When character states for a given OTU were 
underlining. 


APPENDIX II. RAW DATA MATRIX FOR PHENETIC AND CLADISTIC ANALYSES. 
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